i 

ARP  A  ORDER  NO  189-1 

I 


memorandum 

RM-4673-ARPA 

APRIL  1966 


* 


INTERPRETATION  OF 
RECOVERY-TEMPERATURE 
MEASUREMENTS  AT  LOW 
REYNOLDS  NUMBERS 


CLEARINGHOUSE 
FOR  FEDERAL  SCIENTIFIC  AND 
TECHNICAL  INFORMATION 


Hsixloop y 

Uiorofiohe 

%  hot 

S  •S'o 

/; 


pp 


& 


C.  Forbes  Dewey,  Jr. 


n  •  -> 


•^1 


may  2  0 1966 


u  lb 


PREPARED  FOR: 

ADVANCED  RESEARCH  PROJECTS  AGENCY 


K-QflD. 


7 At  K“H  1 1  \J 

SANTA  MONICA  •  C AllYOANIA - 


J 


w  r1**- 


ARPA  ORDER  NO  189  1 


MEMORANDUM 

RM-4673-ARPA 

APRIL  I960 


INTERPRETATION  OF 
RECOVERY-TEMPERATURE 
MEASUREMENTS  AT  LOW 
REYNOLDS  NUMBERS 

C.  Forbes  Dewey,  Jr. 


This  research  is  supported  by  the  Advanced  Research  Projects  Agency  under  Contract 
No.  SD-79.  Any  views  or  conclusions  contained  in  this  Memorandum  should  not  be 
interpreted  as  representing  the  official  opinion  or  policy  of  ARPA. 


■1U 


WWD 


i r oo  «>in  >i  •  iinii  «onic<  >  <  •iiioini* 


DISTRIRITION  STATKMKNT 

Distribution  of  this  document  is  unlimited. 


-iii- 


PREFACL 


An  area  of  great  interest  to  the  ballistic-missile-defense  commu¬ 
nity  is  the  study  of  wake  flows.  Recent  experimental  work  has  been  di¬ 
rected  toward  the  measurement  of  basic  parameters  which  could  be  used 
to  characterize  these  flows.  The  experimental  data  thus  obtained  must 
then  be  compared  with  the  theoretical  data  that  are  available.  The  pur 
pose  of  this  Memorandum  is  to  clarify  the  interpretation  of  the  measure 
ments  so  that  comparison  with  theoretical  results  will  be  more  meaning¬ 
ful. 

This  Memorandum  is  a  part  of  a  study  of  reentry  aerodynamics  for 
the  Advanced  Research  Projects  Agency. 
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SUMMARY 

This  Memorandum  demonstrates  that  the  wire  recovery  temperature 
in  a  steady-state  flow  may  be  accurately  determined  without  a  detailed 
knowledge  of  the  convective  heat- transfer  rate  to  the  wire.  The  anal¬ 
ysis  is  valid  even  for  large  conduction  losses  at  the  wire  supports, 
and  a  similar  analysis  could  be  made  for  transient  conditions. 
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I .  INTRODUCTION 

In  interpreting  recovery- temperature  measurements  in  conventional 

n  2^ 

steady-state  wind  tunnels,  It  is  common  practice'  •  '  to  find  the  aver¬ 
age  temperature  of  a  fine  wire  with  no  electrical  heating  and  to  infer 
the  wire  recovery  temperature  by  a  detailed  accounting  of  the  heat  trans¬ 
ferred  to  the  wire  by  convection  and  of  the  heat  lost  from  the  wire  by 
conduction  to  the  supports  and  by  radiation.  Todisco  and  Pallone^ 
have  recently  proposed  a  variation  of  this  method,  applicable  to  shock 
tunnels,  in  which  the  recovery  temperature  is  determined  in  the  presence 
of  simultaneous  electrical  heating. 

The  purpose  of  this  analysis  is  to  demonstrate  that  the  wire  recov¬ 
ery  temperature  (and  hence  the  local  stagnation  temperature  of  the  stream) 
may  be  accurately  determined  without  a  detailed  knowledge  of  the  convec¬ 
tive  heat-transfer  rate  to  the  wire.  The  analysis  remains  valid  even 
in  the  presence  of  large  conduction  losses  to  the  supports.  This  tech¬ 
nique  may  be  considered  to  be  the  steady-state  analogue  of  the  method 
proposed  in  Ref.  3. 


II.  STEADY- STATE  HEAT  BALANCE 


The  technique  to  be  analyzed  consists  of  two  parts.  Before  flow 
initiation,  the  wire  is  placed  in  a  vacuum  with  a  fixed  current,  and 
the  mean  wire  temperature  is  measured.  Flow  is  then  initiated,  and  the 
change  in  mean  wire  temperature  (mean  wire  resistance)  with  the  same 
current  is  noted.  The  experiment  is  repeated  with  different  currents 
until  no  change  in  wire  resistance  before  and  after  establishment  of 
flow  is  observed.^  The  ideal  recovery  temperature  of  the  wire  is  then 
inferred  from  a  detailed  heat  balance. 

For  simplicity,  assume  that  radiation  is  negligible  and  that  the 
wire  resistivity  varies  linearly  with  temperature  according  to  the  re¬ 
lation 


1  +  cr(T  -  T  ) 

*  - * — —  U) 

R  1  +  or(T  -  T  ) 

where  T  is  a  scaling  temperature  to  be  specified  and  Tf  is  the  temper¬ 
ature  at  which  the  resistivity  coefficient  0/  is  evaluated.  Assume  also 

that  the  wire  thermal  conductivity  k  is  independent  of  temperature, 

w 

that  each  cross  section  of  the  wire  is  at  a  uniform  temperature,  and 
that  the  convective  heat- transfer  coefficient  h  is  constant  across  the 
wire,  i.e. , 


Nu 

0 


hd 

k 

0 


constant 


where  kQ  is  the  thermal  conductivity  of  the  stream  at  the  local  stagna¬ 
tion  temperature  and  d  is  the  wire  diameter. 

The  heat  balance  for  the  wire,  using  the  assumptions  listed  above, 
is 


This  tedious  repetition  may  be  circumvented  in  practice  by  record¬ 
ing  the  wire  resistance  as  a  function  of  current  in  both  a  vacuum  and 
the  flow  and  by  noting  the  intersection  of  the  two  curves. 


-3- 


i2r(x)  +  (^r5)  h  (k„  £ )  ■  nd*h<T  -  tto>  (2) 

where  R(x)  is  the  resistance  along  the  wire  and  T  is  the  recovery  tem¬ 
perature  of  the  flow.  This  equation  is  to  be  solved  with  the  boundary 
conditions  T  =  Tg  (the  temperature  of  the  supports)  at  x  =  ±  (X/2) . 
Define  the  new  variables 


z  »  (x/X) 


s 


_ oT _ 

1  +  o(T  -  T  ) 


★  rs* 

t  =  (T  -  T  )/T 

,  "k  ^ 

where  T  and  T  are  scaling  temperatures.  Then 


R 

=*  1  +  st 

R 


(3) 


(4) 


and  Eq.  (2)  may  be  written 


— £  +  at(l  -  c)  +  a(l  +  bc)/s  =  0 

dz 


(5) 


where 


a 


4i2R*sl  _  4  (i)2  ( 
nd2k  T  n  d'  \  kwX 


(6) 


b 


S(Taw  -  T*)/T  =  *(Tflw  -  T*) /Cl  +  a(T*  -  T.)] 


(7) 
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TiNu  k  i 
o  o 

2 

i  R  a 
r 


(8) 


and  R  is  the  wire  resistance  for  a  uniform  temperature  T  . 
r  r 

In  many  wind-tunnel  investigations,  the  convective  heat- transfer 
rates  are  quite  low,  so  that  the  nondimen sional  parameter  c  is  less 
than  unity.  The  solution  of  Eq.  (5)  is  then  a  cosine  distribution  of 
temperature  between  the  supports,  with  the  wire  temperature  equal  to 
the  support  temperature  at  z  =  ±  1/2. 

*  ~ 

For  convenience,  choose  the  two  scaling  temperatures  T  and  T  and 

the  reference  temperature  T  to  be  equal  to  the  support  temperature  T 

r  s 

(assumed  known  in  a  given  experiment).  Then  s  =  aT  ,  t  =0, 

s  s 

b  =  a^T  -  T  ),  and  the  temperature  distribution  t(z)  across  the  wire 
'  aw  s' 

is 

1  +  s t ( 1  -  c)/(l  +  be)  =  cos  {[a(l  -  c)]^2z}/cos  {[ a (1  -  c)]^2/2} 

(9) 


The  temperature  distribution  without  convection  (c  =  0)  before  flow 

initiation  is  t  (z): 

o 

1  +  stQ  =  cos  (a^2z)/cos  (a^2/2)  (10) 


Integrating  these  results  across  the  wire  to  find  average  wire  tempera¬ 
tures,  t  and  tQ,  yields 

1  +  st(l  -  c)/(l  +  be)  =  2  tan  {[a(l  -  c) ] ^2/2 }/[a(l  -  c)]1^2 

(11) 

1  +  st  =  2  tan  (a1/2/2)/a1/2 

o 


(12) 
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III.  INTERPRETATION  OF  THE  RECOVERY  TEMPERATURE 


The  proposed  experiment  Is  conducted  in  such  a  way  as  to  equate 
the  average  wire  resistances  before  and  after  initiation  of  flow,  i.e., 

K  =  Since  the  resistance  varies  linearly  with  temperature,  this  re¬ 

quires  t  =  t  .  Assuming  that  the  wire  properties  and  the  Nusselt  number 
Nuq  are  known,  Eqs.  (11)  and  (12)  may  be  combined  to  solve  for  the  re¬ 
covery  temperature  Tgw.  After  some  manipulation,  the  recovery  tempera¬ 
ture  is  found  to  be  given  by  the  solution  of  the  following  expression: 


g/(l  4-  be)  =  (tan  3g  -  3g) / (tan  3-3) 


b  =  a(T  -  T  ) 
v  aw  s' 


3  (a/2) 1/2  -  (X/d)(i2R  c*/TTk J)112 

r  w 


g  -  (1  -  c) 


1/2 


(13) 

(14) 


(15) 


(16) 


and  3  and  c  are  defined,  with  R^  and  or  being  determined  from  Eq.  (1) 

and  with  T  =  T  . 

r  s 

The  solution  of  Eq.  (13)  may  be  put  into  dimensionless  form  by  us¬ 
ing  the  measured  mean  wire  temperature  without  flow: 


R 


(17) 


Then  R  *■  R(3,c)  and  in  general  will  depend  on  the  local  Nusselt  number 
Nuq  through  c.  The  parameter  3  is  restricted  to  values  less  than  n/2, 
since  t  is  finite  (see  Eq.  (12)).  The  cosine  temperature  distribution 
has  been  taken  as  0  <  c  <1. 

The  ratio  R(3,c)  has  been  computed  for  pairs  of  values  3  =  0.1, 
0.2,  ...  1.5,  and  c  ■  0.1,  0.2,  ...  0.9.  Somewhat  surprisingly,  the 
value  of  R  is  very  insensitive  to  the  values  of  3  and  c  chosen,  with 
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Rmin  "  1,1995  O  =  0*1,  c  =  0.1)  and  =  1.228  (8  =  1.5,  c  =  0.1). 

An  asymptotic  expansion  for  f32  «  1  yields  R  =  1.2  exactly,  independent 
of  the  value  of  c.  The  interpretation  of  this  result  is  very  important: 
The  recovery  temperature  Tflw  may  be  accurately  determined  without  a  pre¬ 
cise  knowledge  of  the  local  Nusselt  number  Nu  . 

o 

An  a  posteriori  physical  argument  for  the  insensitivity  of  ft  to 

the  local  heat- transfer  rate  may  be  constructed.  The  heat  lost  to  the 

supports  with  and  without  flow  has  been  denoted  by  Q  and  Qq,  respectively 

Without  flow  and  neglecting  radiation,  the  Joulian  heating  i2K  is  equal 

to  Qq.  After  flow  initiation,  the  temperature  distribution  will  become 

more  uniform  across  the  center  of  the  wire,  but  the  temperature  gradient 

at  the  tips  will  increase,  i.e.,  Q  >  Q  .  The  difference  Q  -  Q  must  be 

o  c 

equal  to  the  net  convective  heat  input  in  the  wire  or 

Q  -  Q  ~  Nu  (T  -  T) 
o  o'  aw  ' 

As  the  Nusselt  number  Nuq  becomes  small,  the  difference  Q  -  Q  also  de¬ 
creases  so  that  the  temperature  difference  (T  -  T)  £  (T  -  T  )  re- 

aw  '  '  aw  o 

mains  essentially  constant.  Therefore,  the  ratio  R  is  nearly  independ¬ 
ent  of  Nuq,  and  to  the  accuracy  of  most  experiments,  R  may  be  taken  to 
be  1.2. 

For  this  analysis  to  be  valid,  c  must  be  less  than  unity,  and  a 
steady-state  condition  must  have  been  established.  A  similar  analysis 
may  be  constructed  for  c  >  1  (the  temperature  distribution  follows  a 
cos  h  8z  law)  and  for  transient  conditions  (the  problem  is  linear,  and 
variables  are  separable). 
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